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Abstract 

We show that elastic currents that take into account variations of the tunnel 
transmitivity with voltage and a large ratio of majority to minority spin densi- 
ties of states of the s band, can account for the low voltage current anomalies 
observed in magnet-oxide-magnet junctions. The anomalies can be positive, 
negative or have a mixed form, depending of the position of the Fermi level 
in the s band, in agreement with observations. Magnon contribution is neg- 
ligible small to account for the sharp drop of the magnetoresistance with the 
voltage bias. 
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Tunneling of electrons in metal-insulator-metal junctions is an old phenomenon studied 
from a long time ago . However, it is only quite recently that spin- dependent tunneling 
between two ferromagnetic metals has been shown to produce the magnetoresistance effect 
observed in those systems 0,^. In 3d ferromagnets, most of the spin polarization comes 
from the d bands, while tunneling currents are dominated by s band contributions. This 
is so, because d wave functions are more localized and their effective tunneling barrier is 
higher For Ni, it has been estimated that the tunneling probability of the s electrons 
is of the order of 100 — 1000 times that of the d electrons, thus leading to a positive spin 
polarization in Ni field emission experiments [[7|. In the context of tunneling experiments, 
the large magnetoresistance effect (25-30 %) found in |^,|[ is puzzling, since it points to a 
large polarization of the s band, with a ratio of the densities of states for majority |^A^(*^) (E)^ 

and minority (^N^'^^E)^ electrons at the Fermi level (Ep) of the order of 

Ar(^) (Ef) /Ar("^) (Ef) ^ 2.0 - 2.5 , (1) 

in apparent contradiction with energy band calculations for ferromagnetic metals P]. 

In addition, a remarkable dependence of the junction conductance with the voltage bias 
{V) has been observed at low voltages (of the order of a few hundred millivolts). As usual 
in magnetoresistance experiments, one compares the resistances for the cases where the 
magnetizations at the electrodes are anti-parallel (AP) and parallel (P). In several exper- 
iments reported in Ref. 0,^, the junction resistance drops significantly with the applied 
voltage, with a peak at zero bias (called zero-bias anomaly) that is more pronounced for 
the AP alignment. The effect is also temperature dependent, the peak being less sharp at 
room temperature. Finally, it is found that the junction magnetoresistance (JMR) has a 
large decrease with voltages, up to 60% at 0.5 V in some cases 0. It has been argued 
that this effect can be attributed to the excitation of internal degrees of freedom by hot 
electrons (even at liquid He temperature). Scattering from surface magnons has been pro- 
posed as a mechanism to randomize the tunneling process and open the spin-flip channels 
that leads eventually to a sharp drop of the MR However, this explanation is con- 
troversial, since magnon scattering cross sections are negligibly small to account for such 
a big drop of resistance and no spin-flip events have been observed in experiments with 
polarized injected electrons in tunnelling phenomena [^. Also, the theory given in Ref. 
0] assumes tunneling transmitivities independent of the applied voltages, and uses a per- 
turbation scheme only valid for voltages smaller than ~ 40 mV, while the data extend to 
~ 400 mV. 

In the present Letter, we show that the variations of the conductance with the voltage 
bias can be simply accounted for by the lowering of the barrier height with voltages, as 
given by the Simmons' tunneling theory ^ . The structure at zero bias is obtained, when one 
properly takes into account variations of the density of states with the bias at both magnetic 
electrodes. Assuming that the tunneling current comes from the s band, we formulate a 
simple model with a parabolic dispersion (free-electron like). We obtain different behaviors 
for the zero-bias anomaly, whether the Fermi level is located near the bottom [peak) or top 
of the band (dip). Fitting with the experiments can only be obtained if one assumes a 
large spin polarization corresponding to relation (|ip. 

In order to develope our calculation, one has to rewrite Simmons' formulae with the 
conductance current written in the form 
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J(^)(V) = A^/ dET{E,AsA,V)Nl^\E)Ni,^\E + V)[h{E)-fR{E + V)], (2) 

— OO 

where T (i?, As, 0, V^) is the transmitivity through the barrier for energy E, parametrized 
with the mean barrier height and width As , the index C = P, AP refers to the magnetic 
configuration (parallel or anti-parallel), and Nl^h and fL,R are the densities of states and 
the Fermi distributions for the left and right electrodes, respectively. In ferromagnets, one 
has to distinguish between majority (M) and minority (m) spin bands and the super-indices 
in the densities of states and in the sum in expression (0) label the allowed processes for spin 
conserving tunneling, for both magnetic configuration, P and AP. For parallel alignment, 
the factor of the densities of states that enter in (|) is 

Nt^ {E)N^^^ {E + V)+ iVf ) {E)Nii''^ (E + V) , (3) 

while for the anti-parallel configuration, where majority and minority are interchanged for 
the left and right electrodes, one has to consider 

iV[") {E)Ni,^'^ {E + V) + iVf {E)nP (E + V) . (4) 

Concerning equation (0), several remarks are in order. 

i) In his original treatment of the tunneling problem ||^, Simmons considers the case of 
very fiat conduction bands for the metal electrodes and takes the densities of states 
as constants. However, for s bands the density of states varies as the square root of 
the energy, and for magnetic junctions this cannot be neglected, especially near the 
band edges, where the variation is bigger. Zero-bias anomalies in normal non-magnetic 
metals has been previously reported, in cases where the structure of the density of state 
is important [ pi| . 

a) Expression (|]) involves an integral over all energies, but states that are deep in the 
band are cut off exponentially by the tunneling probability. As a net result, the conduc- 
tance is dominated by electrons that are near the Fermi level, and (0) approximately 
factorizes in the form 

J(^)(V)^| J2 N^''\Ef)N'^^\Ef + V)\ J'^^\V) = D'^^\Ef,V) J'^^\V) (5) 

\a,fj,=m,M j 

where J'^^^iy") is the Simmons' tunneling current as a function of the voltage bias and 

c 

D^^\Ef,V)= N^''\Ef)N^''\Ef + V) . (6) 

a,fi=m,M 

In (^), we are assuming that both electrodes are made from the same ferromagnetic 
metal. The term J^^\V) is the Simmons' contribution, is spin independent and car- 
ries all the information concerning the tunneling barrier. As shown in 0, it has no 
quadratic term in the voltage for small bias, and no zero-bias anomaly. 
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In Fig. 1, we show the variation with voltage of the Simmons resistance for typical 
barriers, with the resistance normalized at zero bias. A large variation is observed in all the 
examples, but the resistance has no peak or dip at zero voltage. Except for the structure at 
zero bias, the overall variation of the Simmons' resistance is of the order of what is observed 
in experiments (or even may vary faster with voltage in some cases). Some experimental 
results are also shown for comparison. 

Next, we introduce the factor D^'-'\Ef, V), defined in (|^), in the conductance calculation. 
We model the density of states of the s bands with a parabolic dependence (free-electron 
like) in the form 



r2 /2me 
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iV('^)(E) = — (^^j ^±{E-E^), a = m,M, 

where Q is the volume of the sample (electrode), rUe is the electron mass, and the ± sign refers 
to the cases where we are in the bottom or in the top of the conduction band, respectively. 
In formulating the Stoner model within a naive band theory, \Em — Em\ should yield the 
exchange of the s band. But Fermi surfaces of transition metals are very intricate, with 
contributions from electron and hole-like carriers and with different shapes for majority 
and minority spin sheets. In this context, Em and Em come from the band structure and 
AE = \Em — Em\ may be very different from the true exchange of the band. 
To parametrized our results, and denoting by Ep the Fermi energy, we define 

rpM — I 771 771 I 

rpm — I 771 77 I 

hip = \1^F — , 

E^' = XE^, A > 1, 

which includes both cases, bottom and top of the band. The ratio of the densities of states at 
the Fermi level is given by n[^^\Ef)/N^i^\Ef) = VX. Several possibilities can be realized, 
wether majority and minority carriers are electrons or holes. When both are electrons or 
holes, the factors D^^\Ef, V) can be expanded in series in V , yielding a linear term in V 
that is responsible for the zero-bias anomaly: 

d'£\v) ^ f[M'")(E^)l'+ [ArW(E^)l') (l±p^P^V\ 



D^^^\V) ^ {2N^"'\Ef)N^^^\Ef)) (1 ± \V\ 



where the ± sign labels the bottom and top cases respectively, with the slopes of the linear 
terms given by 

piP) = 1 



P 



(AP) 



E^ (1 + A) 
A + 1 



AXEf 



When we have a mixed case, i.e. one of the spin is electron-like and the other hole-like, no 
linear term appears in D^^\V). On the other hand, for D'^'^^^V), the slope of the linear 
term is given by 
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where the — (+) sign apphes when the majority carriers are electrons (holes). In Fig. 2, 
we display results of our calculation for examples of typical barriers. The value of the 
magnetoresistance at zero bias was taken from Ref. with 

Nl^'\EF)/Nt\EF) = VX^2.2 . 

In Fig. 2 a), we show the case when the Fermi level is in the bottom of the s band, with 
a linear decrease of the resistance with the voltage bias for both magnetic configurations 
{AP and P). If the Fermi level is in the top of both spin bands, we initially get a linear 
increase of the resistance which, after some voltage value, is dominated by the Simmons' 
term. This case is displayed in part of Fig. 2. In Fig. 2 b), we display the situation 
where the majority band (f) is almost filled (holes) and the minority (J.) is almost empty 
(electrons). The resistance for the P setup, exhibits no linear term. In Fig. 2 a), we also 
show experimental results taken from Ref. We have not tried an optimum fitting with 
experiments, but it is clear that experimental results can only be explained assuming a 
large polarization of the s band. Note that the insets in Fig. 2 a)-2 c) sketch the band 
configurations for both spins. 

The change in tunnel resistance or magnetoresistance (MR) is given by 

Ai? Rap — Rp 
R Rap 

where again, AP and P refer to the magnetic configuration of the ferromagnetic electrodes. 
This ratio, as it is evident from relation (^), is almost independent of the Simmons' term, 
not depending on details of the tunneling process. In Fig. 3 ^4^, we display results of AR/R 
corresponding to the examples of Fig. 2. In B), we take different experimental results found 
in the literature Note that when the Fermi level lies near the top of the band, there is an 
increase of the MR. Eventually, we may reach the minority spin band edge, with a vanishing 
density of states, for which 

Rap ^ oo 

Temperature (T) effects can also be taken into account through relation (^), with the 
broadening of the Fermi distributions, but a rough estimation shows that the effect should 
be similar to that of an applied voltage V ~ 2T, with an effective lowering of the barrier 
height, a smaller resistance, and the softening of the zero-bias anomaly, in agreement with 
experiments. 

From our calculations presented above the following conclusions are pertinent: 

i) The overall variation of the tunnel current with voltage can be explained by elastic 
tunneling using the well known Simmons' formula [Q] and is due to the lowering of the 
barrier by the applied voltage. This is at variance with the calculations in Ref. [^, 
where they argue that this effect is negligible. Therefore, magnons are not needed to 
explain the experiments; 
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ii) The anomalies in the currents and the magnetoresistances can be explained within this 
simple framework, provided that the ratio of majority spin to minority spin electrons 
is of the order of 2.2 — 2.5, for the data of Ref. 0,^]. If one is allowed to choose the 
adequate configuration of the s bands (see Fig.2), a maximum, a minimum or a mix 
of both can appear at the anomaly (as it has been observed in Ref. [|rT|); 

iii) From band structures calculations , it is not clear to us that the above polarization 
of the s band can be justified. There may be other oxidation states inside the metal, at 
the interface, and in the oxide layer, that contribute to the polarization of the current; 

iv) Alternatively, it may also happen, as it has been suggested in Ref. |]6| Jl^ , p!3|| , that the 
current is dominated by conduction paths that provide large values of magnetoresis- 



tance [0 due to domain wall scattering [jl5l, and then there is also contribution of 



d-electrons. In this case, the density of states will have mixed contributions from s 



and (i-electrons, with a variety of topologies in the MR [16 



The main conclusion is that the magnetoresistance is a mapping of the spin up and 
down densities of states in the metals and the barrier and cannot be assigned only to 
the bulk ferromagnetic metals, and many mixing possibilities exist for explaining the 
physical measurements. 
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FIGURE CAPTIONS 



Fig. 1 Variation of the Simmons' resistance with voltage for several tunnel barriers. Data 
is normalized at zero bias. Experimental results from ^ are also shown (solid triangles) as 
a reference. 

Fig. 2 Resistance as a function of the voltage bias for the two configurations of the 
magnetic electrodes and for different s band structures (they are shown in the insets). 
Parameters for the tunneling barriers are given in each figure. Spin | is taken as the majority 
band in all cases. As a reference, experimental results take from [^] are shown in part a), 
where a good agreement with our calculation is obtained. 

Fig. 3 Magnetoresistance, as defined in (|^), for all the cases depicted in Fig. 2. Densities 
of states are adjusted at the zero bias value. In A), we compare with results from while 
part B) compares with Ref. 
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